Calculations of the Local Density of States for Some 

Simple Systems 

Alexandre Losev 

Bulgarian Academy of Sciences, Institute of General and Inorganic Chemistry, 

Sofia 1113, Bulgaria 
and 

Stoyan Vlaev 

Escuela de Fisisca, Universidad Autonoma de Zacatecas, 98068 Zacatecas, Mexico 

February 1, 2008 



Abstract 

A recently proposed convolution technique for the calculation of local densities of 
states (LDOSs) is described more thoroughly and new results of its application 
are presented. For separable systems the exposed method allows to construct 
the LDOS for a higher dimensionality out of lower dimensional parts. Some 
practical and theoretical aspects of this approach are also discussed. 
PACS 71.20.-b, 71.15.-m, 73.20.-r 



1 Introduction 



Recently we proposed and demonstrated a new convolution technique for the 
calculation of the local density of states [1,2]. The idea of this novel approach is 
intuitively accessible, so we had not yet exposed its mathematical description. 
The analytic formulation offers the advantage of a more comprehensive view 
about the interplay of parameters in the utilized models, which could be ben- 
eficial with view to further developments. Exact analytic results being rather 
scarce in this area [3-6], the calculations of the LDOS and other characteristics 
rely on numerical schemes but usually they become a cumbersome task even for 
some of the most simple systems. In consequence, a wealth of alternatives has 
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been produced and among them the ones exploiting tight-binding (TB) mod- 
els appear as particularily useful [7-9]. The TB approximation also provides a 
convenient framework for the present work. The approach considered here pro- 
poses to utilize in a constructive manner the already know LDOSs in analytic 
or numeric form, no matter how they have been arrived at. 

2 Model and Method 

Within the tight-binding model the Hamiltonian of a simple D-dimensional crys- 
tal is written as H = a + 2 PiCos{k.j.^di), i = 1..D. Considered from a statis- 
tician's point of view, the distribution of the values generated by a form such 
as is a 'composition' or 'convolution' [10] of the distributions of its additive 
parts which correspond to the independent variables. Thus, the eigenvalues of 
the system are indeed obtained by a combination of eigenvalues corresponding 
to its components. Accordingly, separating the system into a one-dimensional 
part and a D-1 dimensional subsystem, the LDOS of a the whole system can be 
obtained from the subsystems with smaller dimensions through the recurrence: 

p(D) = p(D-l) * 

where p^"^ = 5 and the sign * denotes the convolution operation. The idea of 
this approach is presented graphically in Fig.l. For a better view all the curves 
in this plot have been scaled to a same maximum height, instead of keeping 
the same (unit) area. If the singularities of the one-dimensional curve are well 
pronounced, they are reproduced in subsequent calculations, as it can be seen 
on this figure. 

In order to have a more substantial derivation we proceed in the opposite 
direction, starting with D = 1. In this case the Hamiltonian is H = a + 
2/3cos(fca:).The distance d is assumed to be d = 1 and the energy e is expressed 
everywhere below through the dimensionless variable E = {e — a)/ 2(3. 

From the meaning of LDOS one can write directly dk ~ p{E).dE i.e. p ~ 
dk/dE. For oscillators, as it has been noted p = dE/dk being the 

group velocity or more generally [11] v = gradk{E). So for this Hamiltonian 
kx = arcos{E) and thus 

the coeiEcient I/tt being written in order to scale the overall density to a 
unit, assuming it is equal to zero outside the energy interval (-1,1). In two 
dimensions H = a + 2(3{cos{kx) + cos{ky)) and in the same way we obtain for 
the partial LDOS 

^ ^l-(E-cos{ky)y 

Integrating for all values of ky and substituting ky with arcos{u), after in- 
verting the integration limits, the LDOS is 

p(2) = fl ^_ du. 

The obtained formula is the canonical expression for a convolution. 
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for which the usual notatfon f = g * h has been adopted. In our particular 

case g and h arc indeed the same function and thus we have a self-convolution, 
which can be of special interest), but this does not need to be so as our examples 
in the next section demonstrate. The same consideration is applicable to higher 
dimensionalities ;is tlic^ multiple integrations are consecutive convolutions e.g. 

would be more generally 

fix) = IJgiix - - vWv)mdridt 

or / = g * h* I for some functions /, g, h and / representing LDOSs. The 
advantages of this approach appear to be fairly obvious: the results for more 
simple systems are utilized to obtain knowledge about more complicated ones. 
Most often such results are available in numeric form and then the convolution is 
implemented simply as the sum = ^ ■ gk_ihi. However a more sophisticated 
procedure might be needed in order to evaluate precisely the integral which 
represents it. In our particular case it may be noted that the quadratures 
for expressions containing 1/-^/ (1 — x^) can be effectively simplified [12]. Also, 
as it is known, a convolution can be expressed in Fourier space by a simple 
multiplication of the transforms and here we have p^^^ = FT{\J^\). Any way, 
the computational advantages of this method can be considerable as it was 
already pointed [1]. 



3 Results and Discussion 

As examples of application we shall consider briefly two cases of semi-infinite 
structures and also offer a comment on the effect of an electric field in the present 
context. 

For the simple model consisting of two infinite wires in interaction or, which 
is the same, an infinite chain of dimers, the convolution method readily offers a 
solution. The LDOS for a dimer is given [13] by the formula: 

p = ^Imag{j^^^), 

where g{E) = 1/{E + ij) is the Green function for an isolated atom and 
t is the interaction between atoms (in the dimer). The explicit expression for 
the LDOS is a sum of two Lorentzians centered at ±t and whose width 7 is 
vanishingly small. In the limit, as they tend to delta functions, the convolution 
produces a sum of two one-dimensional LDOSs centered now at ±t: 

p{E) = \{p^^){E-t) + p^'){E + t)). 

Fig. 2 presents the result for this most simple case. Modifying the value of 
the interaction t would just change the distance between the delta peaks and 
this will be reflected in the overlap of the one-dimensional LDOS, whose width 
is governed by /3. An algebraic treatment for a finite variant of this structure 
is found in Ref.6. More complicated cases are accessible when the poles of 
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the Green function are known. Indeed, if one chooses to consider the LDOSs 
as consisting of Lorentzian (or Gaussian) peaks whose centers are distributed 
following some law, e.g. a cosine for one dimension, then the convolution would 
affect only the distribution, as it produces out of two such peaks a new one with 
a larger full width at half max. 

Another case of a semi-infinite structure is the model of a slab which is 
infinite in along two axes and restricted to only one side of the third. In this 
case the LDOS corresponding to an infinite plane should be convoluted with the 
LDOSs for atoms in a scmi-infinitc chain, which are different according to their 
distance from the origin [13]. This has been illustrated in Fig. 3 which in fact 
reproduces results presented by Haydock and Kelly and who in their turn were 
assessing the possibilities of their own approach by repeating a plot originally 
shown in Kalkstein[14]. The last curve in the plot is essentially a replication of 
the lower row of curves in Fig.l as for atoms far from the beginning of the chain 
the LDOS is the same as in an infinite structure. However one might presume 
that here we have gained some insight into the constitution of these well known 
results. 

Adopting this view, we could reconsider the effect of a constant electric field 
on the LDOS [15]. Here the analytic result is known for the one-dimensional 
case [5] and it has to be convoluted with its an appropriate counterpart in 
order to obtain the results of interest [1]. As it was pointed, a convolution 
can be expressed in Fourier space by a multiplication, and this form makes it 
easier to account for the effect of an electric field: the Fourier Transform of 
Jo{2sin{uj) / F) , where F is the field strength along the axis, matches [16] the 
result for an infinite chain obtained by Davison et al. The FT of the LDOS, 
being represented by a product of such factors corresponding each to an axis, for 
higher field strength along any one of them, turns into a constant (Jo(0) = 1), 
so its original is now a delta function and the dimensionality of the system is 
lowered. The physical aspect of this phenomenon has been considered in the 
work already mentioned as Ref.l5. 

The limitations for this way of proceeding are fairly evident as it relies on 
the additivncss of the Hamiltonians. However some insight in the production 
of the LDOS has been gained. For instance, a p^^' curve now can be seen 
as (the passage to the limit for) a weighted sum of 'U' shapes centered at 
different energies. If second neighbours are to be included in the tight-binding 
Hamiltonian, it ceases to be additively separable but it still can be factorized, 
now as a product, and then the summation would include not just weighted and 
displaced one-dimensional 'U' shapes but ones with a different spread. This 
suggests that our constructive approach could perhaps offer further possibilities 
for the calculation of more complicated cases. 
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Conclusion 



The proposed technique allows to utilize effectively already obtained results 
and to extend them further. It is able to generate new ones and in many in- 
stances otherwise difHcult to obtain results became easily accessible. However 
the method relies on the separability of the Hamiltonian, which severely limits 
its scope. The analytic formulation of the method suggests a novel view even if 
the complexity of the calculations most often defies their completion. The con- 
volutional form has a rather intuitive meaning, which might be more intelligible 
than the abstract elliptic integrals expressing the same result. 
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Figure Captions 

Figure 1. The LDOS of a 3D simple cubic crystal obtained by self convolution. 
Upper panel: the convolution of two ID LDOS produces a 2D LDOS; Lower 
panel: the convolution of the 2D LDOS and ID LDOS produces the 3D LDOS. 
Figure 2. LDOS for an infinite chain of dimers obtained by convolution. 
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Figure 3. LDOS in a semi-infinite slab for the first three atomic layers and 
bulk (z = 0, l,2..oo). The dotted curves present on a smaller scale the LDOS 
in a scmi-infinitc chain while the inset shows the 2D LDOS in a plane used to 
produce the curves. 
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